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Abstract 

We study the time evolution of small classical perturbations in a 
gauge invariant way for a complex scalar field in the early zero curva- 
ture Friedmann-Lemaitre universe. We, thus, generalize the analysis 
which has been done so far for a real scalar field. We give also a 
derivation of the Jeans wavenumber in the Newtonian regime starting 
from the general relativistic equations, avoiding the so-called Jeans 
swindle. 

During the inflationary phase, whose length depends on the value 
of the bosonic charge, the behavior of the perturbations turns out to 
be the same as for a real scalar field. In the oscillatory phase the 
time evolution of the perturbations can be determined analytically as 
long as the bosonic charge of the corresponding background solution 
is sufficiently large. This is not possible for the real scalar field, since 
the corresponding bosonic charge vanishes. 
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1 Introduction 



The developments in particle physics and cosmology suggest that scalar fields 
may have played an important role in the evolution of the early universe, for 
instance in primordial phase transitions, and that they may constitute part of 
the dark matter. Moreover, scalar fields are an important ingredient in most 
of the particle physics models based on the unification of the fundamental 
forces, as for instance in superstring theories. Scalar particles are needed 
in cosmological models based on inflation, whose relevance is supported by 
the results of the COBE-DMR measurements that are consistent with an 
Harrison-Zel'dovich spectrum |IJ. 

These facts, in particular inflation, motivated the study of the coupled 
Einstein-scalar field equations to determine the time evolution and also the 
gravitational equilibrium configurations of scalar fields, which may form so- 
called boson stars 0. [|. 

A detailed study of the solutions of the Einstein equations for a homo- 
geneous isotropic Friedmann-Lemaitre universe with a real scalar field has 
been carried out in particular by Belinsky et al. M, 0] , Piran and Williams 
0. The corresponding thorough investigation for a complex scalar field has 
been performed in Ref. M 

The knowledge of the behavior of the background solution in the early 
universe is of importance for the analysis of the time evolution of the per- 
turbations [0, TO]. This is of relevance if scalar fields make up part of the 
dark matter and form compact objects, such as boson stars, or trigger the 
formation of the observed large scale structures in the universe. 

In this paper we study the time evolution of small perturbations in a 
gauge invariant way for a complex scalar field in the early zero curvature 
Friedmann-Lemaitre universe thus, generalizing the previous analysis for a 
real scalar field We consider a potential with a positive mass and a 



quartic self-interaction term and closely follow the treatment developed by 
Mukhanov et al. [ID|. During the inflationary phase the behavior of the 
perturbations is similar to the one of the real scalar field. This is not sur- 
prising, since along the separatrices the phase of the complex scalar field 
remains constant and thus inflation is essentially driven by one component 
of the field ||. The short wavelength perturbations get smeared out, whereas 
the long wavelength perturbations increase. On the other hand, the metric 
long wavelength perturbations are entirely determined by the scale factor 
and hence survive the decay of the complex scalar field background solution 
after inflation. 

We extend our investigation of the behavior of the perturbations also to 
the oscillatory phase, which takes place after inflation. Contrary to the real 
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scalar field for which the bosonic charge of the background solution van- 
ishes, the charge is different from zero in our case. If it is sufficiently large, 
we can even find analytically the time evolution of the perturbations for a 
massive complex scalar field. It turns out, that the long wavelength pertur- 
bations decrease with time, whereas the short wavelength perturbations, at 
best, oscillate around a constant value. Therefore, we do not expect the for- 
mation of gravitational seeds around which structure can form for the case 
of a massive complex scalar field with a sufficiently large bosonic charge. 
The other cases have to be treated numerically. In order to put constraints 
on the parameters of the scalar field potential from the measured cosmic 
microwave background (CMB) anisotropics [|lj, we have to analyze the quan- 
tized perturbations. Here, we have restricted the analysis only to the classical 
perturbations. 

Moreover, we also derive the Jeans wavenumber in the Newtonian regime 
starting from the general relativistic equations, thus avoiding the so-called 
"Jeans swindle" . For the case where there is a quartic self-interaction term 



in the potential the result we find differs from the one given in Ref. [12 . 

The paper is organized as follows. In section ^| we present the basic 
equations, which we will use. In section [3] we give the derivation of the Jeans 
wavenumber starting from the general relativistic equations. Section |] is 
devoted to the analysis of the time evolution of the classical scalar mode 
fluctuations in a fully general relativistic way, both in the inflationary and in 
the oscillatory phases. The numerical results are presented in section [5| and 
in section |] we present briefly our conclusions. 



2 Basic equations 

We will consider the linear scalar mode perturbations since they are the only 
ones which contribute to the energy density fluctuations [JIB] . As usual, we 
expand the scalar perturbations in terms of a complete set of harmonic func- 
tions Yk, which are the eigenfunctions with eigenvalue — k 2 of the Laplace- 
Beltrami operator A defined on constant time slices S^. In the following, we 
will omit for simplicity to write the subscript k from Y. 

The zero curvature Friedmann-Lemaitre metric — including the first order 
perturbations — is then given by (T^j 

d s 2 = 9flv e^e v = -e°e° + <%w, (i) 
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where 



and 



0° = adr] = a(l + AY) drj , 
= (?dri + a[(l + H L Y) dx i + H T Yj dx j 
= -aBY i dri + a [(1 + H L Y) dx i + H T YJ dx J 

Y% = — k Y|j , Y^ = k Y\ij + — Y5ij . 



a is the lapse function and (3 l are the components of the shift vector. A, 
B, Hl and Ht are, respectively, the amplitude of the perturbation of the 
lapse function, the shift vector, the unit spatial volume and the anisotropic 
distortion of £,,. All these functions, as well as the scale factor a, depend only 
on the conformal time rj. The subscript u denotes the covariant derivative 
on (E^, 7), where a 2 7 is defined as the unperturbed induced metric on E^. 
It is easy to verify that Y^ is traceless. 
We take for the matter action 

S m = f( [J<r (e M $ e,$* + e u $e^*) + V (<t>, $*)] , (2) 

where ( is the volume form and e M is the dual basis of 9^. We expand the 
complex scalar field $ to first order, where 0(77) is the zero-order background 
solution and 5(j)(r})Y is the first order term. The potential V is such that 
S m is invariant under a global U(l) symmetry. Thus, V has to fulfill the 
following conditions 

dV 



dV 



dV 



,d 2 V 



*2 



d 2 V 



, d 2 v 

d 2 V 



= 0, 
= 0, 

= 0. 



d<p* d<p* 2 d<pd<p* 
The C/(l)-global symmetry of S m implies a conserved current 



where 



J(0) = 

5J = 



a 
ik 



eo , 



+ 



-J^A\Ye 



[050* - 0*50] Y*ei . 



(3) 
(4) 
(5) 

(6) 
(7) 

(8) 
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Dot means the derivative with respect to r\. By varying S m with respect 
to the metric, we obtain the energy-momentum tensor up to the first order 
terms 



T = ee <g> e + e <g> S + S ® e + T = T (0) + 5T, 



(9) 



where 



e<°> + Se 
1 



+ 



2AY<j>*4> + a 2 ^50 + a 2 ^-5<p* 



(10) 



s 

T 



S(°) + 5S = — (<p5<f)* + <p*5<f)) Y i e l 
T (o) + ST 



(11) 



V - a 2 V + 



,dV 



dV 



— a 



+ 



\Y 



;i2) 



The background quantities and the first order perturbation components can 
be easily identified. To derive the Einstein equation using the 3 + 1 formalism 
[ |14| , we first have to compute the second fundamental form K for the slices 
Ejj, for which we get 



K 



1 



hS ij + ka Y ij + hn a Y5 l] 



6; ® e 



3 > 



(13) 



where h = | is the Hubble expansion rate 



On 



-B 



and 



, kB H L 
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The Einstein equations are then given by |TJ 

R + [Tr(K)] 2 - Tr(K 2 ) = ISiiGe , 
V ■ K - VTrK = 8ttGS , 

afRic(g) + 2K 2 



(14) 
(15) 



(16) 



d t K =LpK - Hess(a 

+ KTrK - 8ttGT - AnGg(e - TrT) 

where Tr denotes the trace, Hess the Hessian, g the induced metric on £ 
R and Ric(g) are, respectively, the Ricci scalar and the Ricci tensor of the 
slices Y> v . The background solution trivially satisfies eq.(15), whereas from 
eq.(14) and eq.(16) we obtain 



v 



8ttG 



* + a 2 V 



3h 2 

h + 2h 2 = 8nGa 2 V. 
^From eq.(16) we get for the terms in first order 



(17) 
(18) 



fc%-($ ff + V) =-Y8i. 



—ip — 5hip k 2 ip 



I k 2 
+ h$ g + ( 16irGa 2 V - — 



8iiGa 2 \^-5ip 



(19) 



dV 



6<p* 



where we have introduced the gauge invariant quantities ]T6|, [13], [IU| 



if, = 


H L + ~H T - -ha 9 , 


(20) 




50-0^, 


(21) 


% = 




(22) 






(23) 



For % ^ j we obtain 



The meaning of the gauge invariant quantities is obvious on small scales 
(i.e. k ^> 1), in which case they reduce to Hl + \Ht (the intrinsic scalar 
curvature perturbation), <50 (the complex scalar field perturbation) and A 
(the gravitational potential), respectively. 
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The first order terms of eqs.(14)-(15) combined with the trace of 
eg . (|19"D lead to the following Einstein equations 



6hip + (2k 2 + 6h 2 ) i\) = 8tcG U>*6(p + jri'v* + 2 #*^ 

< 24 > 

^ + hip = -4ttG U*5ip + <j)5ip*) , (25) 

— 6<p + 7-— 6<p* ) , (26) 

0(p d(p* J 

where we used eq. (p3|) . 

By varying S m with respect to $* we obtain the Klein-Gordon equation, 
which to the zero order is given by 

dV 

<P + ^<P + a 2 — = 0. (27) 

By directly expressing the first order perturbed Klein-Gordon equation in 
terms of the above defined gauge invariant quantities, we get 

... ■ ■ dV d 2 V d 2 V 

5<p + 2h5<p + k 2 5ip + 40V - 2a 2 — ^ + a 2 -—5^>* + a 2 — —5<p = 0. 

ocp* d(p* z d(p*0(p /2g\ 

Similarly by varying S m with respect to $ we obtain the complex conjugate 



of eqs.(27)- (p8|) . It should be noticed that the first order terms are entirely 



determined by eqs.(p4|)-(p5|), eq.(p8|) and its complex conjugate. These equa- 
tions hold for all values of k and, therefore, the range of validity is not limited 
by the horizon scale. The solutions of these equations will be valid as long 
as \5<p/(f>\ and ip are sufficiently small, irrespective of k. 

The background solution is also established by eq. (O) , eq.(P7p and its 
complex conjugate. The fact that eq. (|l8|) and eq.(|26f) are automatically sat- 
isfied is due to the Bianchi identities. 

The conservation of the current J leads to a set of two additional equa- 
tions. For the background part of the solution, we get the conservation with 
respect to the conformal time of the bosonic charge 



„ i a 2 



(29) 



2 

and for the first order term the expression 

d v £ = a 2 k 2 (<f>*5(p - <t>5<p*) , (30) 
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where 

f = 4a 2 ^ (0*0 - 00*) + a 2 (0*5^ - <j>5ip* + 0oy?* - <f>*5<p) . (31) 

For very large wavelengths £ is also a conserved quantity, since then the right 
hand side of eq.(|30D vanishes up to first order. 



3 The Newtonian regime 



In this section, we derive the Jeans wavenumber starting from the general 
relativistic equations. This way, we avoid the so-called "Jeans swindle" p~7 . 



The gauge invariant formalism is very well adapted to this problem, since in 
the longitudinal mode {B — E — 0) ip is just the Newtonian potential. 

To obtain the Jeans wavenumber we have to derive the dispersion relation. 
We consider a potential with a mass and a quartic self-interaction term. The 
Newtonian regime is realized when the expansion of the universe can be 
neglected and the horizon size is much larger than the wavelength of the 
perturbation (i.e. k ^> h and h ~ e, where e<l). With these conditions it 
follows from eq. (|i~7l) that 



k 2 >87rG00* 



and 



k 2 >87rGaV = 87iGa 2 m 2 00* + A(00 



£ *\2 



(32) 



(33) 



We obtain two conditions from eq . (|17D , because the right hand side of this 
equation is a sum of two positive terms. ^From eq.(|T8"D we have h k 2 and 
thus I a/a I k 2 . We have now to find the solution for the background in the 
Newtonian regime. Inserting 



into eq. 



Imposing 



we get 



4> = f/a 



f + f a 2 m 2 + 2Xff* -h-h' 



0. 



a 2 m 2 + 2A//*> h + h 2 =\a/a\ 



(34) 
(35) 
(36) 



— which means, that for an oscillation period of 0, a can be considered as a 
constant — eq . (|3~5D can be solved and we get 



/ = foexp(iur]) 



(37) 
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where 



\la?m 2 + 2A/ 2 , 
and /o is an integration constant. Using eq. (|34|) , we have 

4> = b exp(iurj) 

with b Q = fo/a, which within our approximations is constant. Hence, 
eqs.(32)-(33) and eq. (|36|) can be rewritten as 



(3? 



(39) 



k 2 > 8nGb 2 a 2 (m 2 + 2\b 2 
k 2 > 8irGb 2 Q a 2 (m 2 + Xb\ 



< a 2 m 2 + 2\a 2 bl. 



(40) 
(41) 

(42) 



To obtain the dispersion relation we perform the following expansion 

Y 



hY\ (. .tpY\ 

b H exp tour] + % 

a \ a J 



(j)+-(b 1 + ib Q ip) e 8 ^ + 0(2) , 



/ * 

v> = — . 

a 



(43) 
(44) 



As next, we rewrite eq.(|2"4]) and eq.(^) in terms of &i and Eq.(|2"8]) will 
then split into a real and an imaginary part. With the following expansion 



&i = b 10 exp (Q77) , ip = v?i exp (^77) 



^ = ^10 exp (fir)) 



(45) 



and using the above mentioned Newtonian approximations, we get a set of 3 
equations 

(46) 



(47) 
(48) 



l6TrGa\b 10 [m 2 + 2Xb 2 ) + ^Gb 2 Q VLuj^ w - k 2 ^ 10 = 
610 (ft 2 + k 2 + 4Aa 2 6 2 l ) - 2Qubo(p 10 - 2^i a 2 & (™ 2 + 2Xb 2 



+ fc 2 ) + 4*io6 ^ft = . 
After some algebraic manipulations, we find 



167rGb (jj 2 





(4uj 2 + n 2 + k 2 + 4\a 2 b 2 ) a 22 
a 33 



8irGb 2 cutt 



-k 2 





'b w ' 








*10. 



(49) 
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with 



d22 



-327rGb u 4 + k 2 Q 2 + k 4 + AXa 2 b 2 k 2 



(50) 



and 



1 



«33 



-MnGb 2 u 4 k 2 + 2k b + 8Xa%k 



21.21,4 



327rGb a 2 uj 2 k 2 
+ Q 2 (2k 4 + 64vrG6^ 2 + 2k 2 ) Uu 2 + Q 2 + k 2 + 4Aa 2 6 2 l 



(51) 



033 = is the dispersion relation, which is the required condition for having 
a non-trivial solution of eq.(49). The Jeans wavenumber, corresponding to 
the solution of the dispersion relation with Q = 0, is given by 



k, = 2a 2 bn 



X 2 b 2 + 87rG(m 2 + 2\b 2 ) 2 - A6 C 



(52) 



For A = we obtain the same result as in Ref. 12 . Notice that, since we 



have taken the conformal time, we have an over-all a 2 factor. For A 7^ 0, we 
disagree with the expression for the Jeans wavenumber givenQ by eq.(31) in 

Ref-H 



This is due to the fact that eq.(l) in Ref. [0 is not compatible with 
a quartic self-interaction term. Indeed, for the mass-less complex scalar field 
with A 7^ 0, we expect to have a non- vanishing Jeans wavenumber as it is the 
case for our result, given by eq.(^), but not for the corresponding eq.(31) 
in Ref. [0 . Only disturbances whose wavenumber is smaller than the Jeans 



wavenumber kj can grow. Hence, the knowledge of the Jeans wavelength 
Xj = 2njkj gives a rough idea of the size of the objects which can be formed 
by gravitational collapse. What we do not know, however, is the rate at 
which the perturbations will grow. 



4 Classical perturbations of a complex scalar 
field 

In this section we analyze the time evolution of the scalar mode fluctuations 
in a fully general relativistic way. We will consider the long wavelength and 
the short wavelength limits separately. To solve the Einstein eqs.(^4|)-(^), we 
first define the complex valued function U (rj) as the solution of the following 
differential equation 

U + hU = -47iG<p*5<p . (53) 

1 7,2 8TTGm 2 (m 2 + \\\ 2 a 2 ) 
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For \5<p\ < \(f)*\/4irG the complex scalar field perturbation is solution of the 
constraint eqs.(|24])- (|25|) if and only if the Bardeen potential i[> and U satisfy 
the following equations 

$ = U + U* + ^ , (54) 
a 



U + 2 ( h - t- ) U + ( k 2 + 2h - 2hC) U = (A - // - ir 
\ 6* \ 6* I 2a 



(55) 



the r.h.s of the last equation represents a source term and k^ is a real in- 
tegration constant. In fact the term k^/a in eq.(54) is a solution of the 
corresponding homogeneous eq. (|25|) , namely ip + hip = 0. In the real scalar 
field case there is no source term fllP] . To simplify the notation we omit 
to write the explicit k dependence of U. Setting U = — , eq.(55) can be 
rewritten as 



u + k 2 u + (h-h 2 + ^-2^)u = -^(k 2 + h-h 2 ) . (56) 



1*2 

^ i, 

20 

With g = -j; and using the background equation we find 



u + k 2 u+(^^-^)u = -^(k 2 + h-h 2 ) . (57) 
V a 2 gj 20* V ! 

By solving this last equation we get immediately U. Using eq.(54) and eq.(|53|) 
we obtain if> and Sip, respectively. 

In the analysis of eq.(|5TD we will restrict ourselves to long and short 
wavelength perturbations in the inflationary and the oscillatory phase of the 
background solution. 

4.1 Perturbations during inflation 

We consider inflation generated by the potential 

V = m 2 (jxj)* + A(00*) 2 , (58) 

for which during inflation |0| 2 $ is asymptotically zero ||, where we have set 
(f) = |0|e i1? . ^From eq.(^) we get 

a 2 |0| 2 tf = H, (59) 

where S is the constant bosonic charge. Hence, inflation must start with a 
small value of E/a 2 . On the separatrix, where inflation occurs, E/a 2 decays 
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exponentially, whatever the values of m and A are. As a consequence eg. flSTp 
reduces to 



i 7 2 ^ 

M + K U U 

9 



20 



t (k 2 + h-ti 

A* V 



(60) 



Moreover, we see from eq.(^) that for small values of H the inflationary 
phase lasts longer. Eq.flSTjD is actually valid not only for the inflationary 
stage but also when E/a 2 is negligible and |0*|>47rG \5(p\. 
During inflation we have for A ^ and A = 0, respectively 



a 2 ^ 



and 



g/g 



-a 2 V, 



-a 2 H t 



-2a 2 H t 



2 h-h' 



where we have denoted 



9V 



(61) 
(62) 

by and introduced H = ~ = ^ with t 
being the real time defined as c?t = a CZ77 (i^ means derivative of H with re- 
spect to real time). This result can be easily obtained using the inflationary 
background solutions (see ||). Since inflation occurs only when V ^> 4>t4>ti 
we have \H t \ <C H 2 . This means that g/g <C a 2 H 2 and g/g increases expo- 
nentially. Hence the time interval, where k ~ g/g, is very short. 

For k g/g — the long wavelength perturbations — the integration of 
eq.flSOD gives 



u 



h9+k ^Sv 2 - g Sj 2 1 



rv kjp (h 



20* 



drj 



drj 



(63) 



with k\ and k 2 being complex integration constants. Using the background 
equations and the fact that the phase of the complex scalar field is constant, 
we can rewrite eq.(|63|) as 



koe 



-2i-d 



d 



k v 



u = -f- - / a 2 drl - . (64) 

8ttG0* drj [a J J 20* 

In this last expression the second integration constant has been absorbed in 
the integral. It follows that for this case the second term in eq.(64[), which 
is due to the source term (see eq.(55)), gives no contribution as we will see 
below. With eq.(B3) we obtain for the gauge invariant metric perturbation 



4, 



-Re(k 2 ) 
2Re (k 2 ) 



a 2 drj 



2Re ( h 



d_ 
~di 



a dt 



H , 

1 / adt 

a 



2Re k 



H 



-1 



H 



-1 



H 



-1 



H 



-1 



H 



-1 



H 



-1 



(65) 
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and during the inflationary stage 



ijj ~ -2Re (k) 



El 
if 2 ' 



(66) 



where k 2 = k 2 e 2t ® /8tcG and Re means the real part. For the gauge invariant 
field perturbation we get 



5 ip ■ 

and during inflation 

5ip 



k 2 d>* r> 



AnGa 2 



a 2 dfj 



k^t ft 



AnGa 



a dt 



AnGH 



(67) 



(6* 



For k ^> g/g — the short wavelength perturbations — eq.(|60|) can be 
integrated and we obtain 



u = 71 cos krj + 72 sin kr] + 



cos kr] fi k^k 2 



k 



sin krj drj 



sin kr] p k^k 2 
k 



(69) 



cos kr] dfj 



As a consequence, 



i/j = 2Re 
~ 2i?e 



(71 cos(kr]) + 72 sin^r?)) 



a 



3k 2 



4M£\ + 3k 2 j 



t (71 cos(k J 72 sin f A; ^ 



a 



(70) 



1 f fc 

^ = ~A^G\ a COS ^ 77 ^ ~ 7l * sin (^)] 

( (71) 

H ^- [7* cos(kr]) + 72 sin(A;77)] 

where 71,72 are complex integration constants and M 2 = l/8nG. Notice, 
that also for this case the term proportional to k^ does practically not affect 
ip. During inflation for a perturbation wavelength smaller than the Hubble 
radius — k 2 > a 2 H 2 ~ 0*/20* — we find 

k 

6if ~ ~A~^G~a ^ 2 cos ( kr] "> ~ 7 * sin (^)] • ( 72 ) 
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The slow-rolling approximation, required for having a sufficiently long in- 
flationary stage, leads to a slowly variation of the mean value of ip. The 
behavior of 5ip is governed by the 1/a factor, which decreases rapidly. It fol- 
lows, as expected, that the short wavelength fluctuations of the scalar field 
are smeared out. 

We consider an initial perturbation with wavelength smaller than the 
Hubble radius (k 2 ^> a 2 H 2 ^> g/g), which will be outside it at the end of 
inflation (k 2 <C a 2 H 2 ). Indeed there are wavelengths that fulfill these condi- 
tions, since a 2 H 2 increases exponentially (see Fig. 1). At the end of inflation, 
the evolution of the gauge invariant metric potential is given by eq.(65). Later 
on, when the universe is dominated by relativistic particles, the scale factor 
is a power-law function a oc t^. Hence, the Hubble radius increases more 
rapidly than the fixed comoving wavelength. The metric perturbation can 
thus re-enter inside the Hubble radius and induce fluctuations on the ordi- 
nary matter. These fluctuations would be the source for the anisotropics in 
the CMB, which have been measured by COBE [Q. 

At the Hubble radius crossing, with eq.(65), the metric perturbation is 
given by 

i> = 2Re (k 2 ) (73) 

We can now use the relation between k and Sip at the time t c , defined such 
that k = g/g, and substitute it into the last equation. We get 

$ = - Re ( 6 -^He~ 2i A -J— . (74) 

There is an amplification of the metric perturbation between the time it leaves 
the Hubble radius and the time it re-enters. Using the asymptotic solutions 
for the inflationary stage this amplification can be derived analytically. As 
an example for A = we obtain 

ip in _ m 2 t 2 c 



4>out 3 (ji + 1) 

This generalizes the result of previous investigations for the real scalar field 
case [pll| . In eq . (|74]) , the gravitational potential still depends on the value of 
the complex scalar field perturbation at the time t = t c . In order to be able 
to put constraints on the parameters of the model, by comparing with the 
observed CMB anisotropies, one should be able to express ip as a function 
of the background quantities only. To achieve this one has, as for the real 
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scalar case |18[], to quantize the complex scalar field. However, this task is 
not as simple as for the real scalar field, since in the lagrangian there is an 
additional interaction term between the real and the complex part of the 
scalar field. We plan to come back on this issue in a future publication. 



4.2 Perturbations during the oscillatory phase 



After the inflationary stage the background solution goes through an oscil- 
latory phase. This corresponds in the phase portrait to winding around a 
focus point A for which = = 0. For a massive complex scalar field 
the asymptotic behavior of the background solutions around the point A is 
discussed in Ref. ||. Expressing the asymptotic behavior as function of the 
conformal time, we have 



3a m 2 r/ 2 



(76) 



2M„ 



v^3 a m 3 r] 3 



cos $30 cos(a m 3 ?7 3 — 771) 
+i sin $30 cos(aom 3 r] 3 — 772) 



(77) 



2^3 M„ 



cos $30 sin(a m 3 ?7 3 — r/x) 
+i sin $30 sin(aom 3 ?7 3 — r] 2 ) 



(78) 



for 77 — > 00. ao, T]i, r\ii $30 are dimensionless integration constants. Due to 
the following scaling behavior 

«0o lfi v) = al /3 a(v) , 

H a o 1/3 v) = <f>(v) 



and 

0(a~ 1/3 r/) = a o /3 0(r/) 

of the background equations, we can restrict the analysis to ao = 1. 

For the study of the first order fluctuations eq. (|56|) can apriori only be 
used as long as |0*|>47rCr \5<p\. |0*| is a damped oscillating function which 
has its minimum on an oscillation period when 

b = cos 2 $30 sin 2 (m 3 r/ 3 — rji) + sin 2 $ 30 sin 2 (m 3 ?7 3 — r/ 2 ) (79) 
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is minimum. If b is large enough (such as b/rj > \S(p\) then eq.([5lj) will be 
valid for a large number of oscillations. On the contrary, if the minimum 
of b is very small, then there will be for each oscillation period two time 
intervals where eq.(^) is no longer expected to be valid. One can easily see 
that the minimum of b vanishes for $30 = mr/2 or (r] 2 — rji) = nir, where n 
is an integer. These conditions correspond to a vanishing bosonic charge S. 
Indeed, inserting eqs.(|76|)-([78|) in eq.([29l), we find 



cos $30 sin $30 sin(ry 2 - r]i) 



10) 



where S max = 18M 2 m is the maximum bosonic charge. The maximum value 
of all possible minima of b is reached for $30 = 7r/4 + mr/2 and 772 — ^71 = 
7r/2 + Itt (/ being an integer). For these and only these values of $30 and 
{V2 — Vi)> b = 1/2 and |H| is maximum. One can show that by increasing 
the minimum of b also the absolute value of the bosonic charge increases. 
Following these lines, we see that the time interval, where eq. (|56D is valid, 
increases with |S|. Hence, we will restrict our analysis to 5 



With eqs.([76| ) - ([78|) we can easily show that 



9m 6 r) 4 



(b-2f 



^2 



772 

' — 'max , 



The right hand side can never be zero for any value of b. It follows that 



Thus, eq.(j56|) reduces to 



\h-h z \ < 



P 



u + k u—'—u 
P 



%(e + h- h 2 ) 



(82) 



33) 



where p = 1/0*. There are three regimes to study. The first is given by the 
condition k 2 3> \p/p\ 3> \h — h 2 \, the second by \p/p\ ^> k 2 ^> \h — h 2 \ and 
the third by \p/p\ 3> \h — h 2 \ 3> k 2 . 

For the short wavelength perturbations — k 2 ^> \p/p\ - eq. (|69D can 
again be used. Using that for 5 ~ H ma:r 



\/2Mp ^±i(rn 3 r) 3 -r]o) 

\f?>rrfirf' 



and 



VSM p ±i( r n 3 rf i -Vo) 



35) 
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we get that 



4> 



2Re 
2Re 



(71 cos(krj) + 72 sm(kr])) 



+ + 0(1/^ 
a 



t ( 7! cos ( k J — J +72 sin I A; J 



a 



+ 



(86) 



and Sp is given by eq.(71), since it can be shown that the inhomogeneous 
solution is of order 1/rf 1 . The different signs in eqs.(84)-(85) depend on the 
integers n and I. We see that the dominant term of ip decreases like 1/a. 
With eq.(71) we have 



8p> 



AttG, 



[7J cos(krj) + 72 sm(kr})] 



57) 



since in general 



CL(f)* 



m r 
2b 



cos ^30 sin (2m rj — 2r\ 



+ sin ^30 sin(2m 77 — 2r\<i) — i— 



(88) 



and fc/a is a decreasing function. It follows that the real and the imaginary 
part of 5ip oscillate around a constant value. 

For the long wavelength case — k 2 \h — h 2 \ <C \p/p\ — u is now given 
by eq.(|63|), where we have to substitute g by p. It is easy to see that for 
large t 



t 2 



where 



k-i = 2k 2 M 2 (sin 2 $30 — cos 2 $30 + 2% sin $30 cos 



'30 



59) 



(90) 



It should be noticed that for E A; 3 + fcg is a very small real number. 

In the long wavelength regime the metric perturbation is essentially driven 
by the expansion. For the computation of Sip we use also eq. (|53|) and obtain 



Sip 



k, 



AttGcl 



AttG 



(91) 



For this case the source term gives only a negligible contribution to 8<p, and ip 
gets another term proportional to 1/a, which is absorbed in its homogeneous 
solution. 
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In the intermediate case, where \h — h 2 \ <C k 2 <C \p/p\, u is given by 



u = kxp + k 2 p — -v \ ^ \ ~r~ dr l 
j pz j pz j 2(f)* 



<"n k^k 
20 



d?7 . (92) 



After some computations we see that the solution of ip is also given by eq. i 
For the complex scalar field perturbation instead we have 

Hence, the first order perturbation of the scalar field decays in the long 
wavelength and in the intermediate regime, whereas it oscillates around a 
constant value in the short wavelength. Therefore, during the oscillatory 
phase there can be no formation of gravitational seeds, like boson stars, by 
a massive scalar field with sufficiently large bosonic charge. 

Using the background asymptotic solution ||, we obtain the analytic 
behavior of the metric and the complex perturbations. They are reported in 
Table 1. For all the other cases the solutions have to be found numerically. 



5 Numerical results 

To complete our analysis, we study the evolution of the perturbations after 
inflation, which implies that the bosonic charge S is close to zero. 

As a first step we integrate the background equations from an initial 
singularity point S, where the equation of state is given by = (stiff 
matter) and the asymptotic behaviors of 0, and H are as follows M 



M p V3 



C+i-\n(t/t ) 



(94) 



M p V3 



CmH i 
2~ + 3t 



(95) 



H = — 

3t 



(96) 



for t — > + , C and to being integration constants. It should be noticed that 
the bosonic charge S depends directly on the choice of the value of C and 
to- Hence, starting from S we can freely fix the value of H. In order to 
have a long inflationary phase we must choose a small value of S, which is 
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achieved for instance by setting C = 1 and to — 100. With these values 
the background starts with a stiff matter regime and, after an inflationary 
phase {h exponentially increases, see Fig. 1), ends in the oscillatory phase 
(h oscillates around a slowly decreasing curve, see Fig. 1). For a massive 
complex scalar field inflation is also less effective when A increases (indeed, 
the maximum of h decreases with A, see Fig. 1). This, because the potential 
becomes less flat and as a consequence the slow-roll approximation, needed 
for inflation to occur, is no longer well satisfied. 

As a next step we integrate numerically eq . (|26|) , eq. (|28"D and its complex 
conjugate using the background solution. We control our numerical compu- 
tation by looking if eqs.(|24])-(|25D propagate. Moreover, we check if we recover 
the analytic solutions described in section f|. 

The Bardeen potential ip behaves as predicted during inflation and re- 
mains, up to a smooth oscillation, constant after inflation. This behavior is 
valid for any value of the parameters k, m and A. 

For a long period during inflation \6<p\ decreases and thus, at first glance, 
does not behave as expected from our analytic analysis (see Fig. 3). For 
the cases A = 0.1 and m = with A = 1 we even get that during inflation 
\8<p\ is only decreasing. This may seem to be in contradiction with the an- 



alytic results obtained in section ^J]. However, this is not true since we do 
not start the integration of the background right away from the singularity 
point, where inflation occurs. If we did so, we would have obtained the ex- 
pected analytic behavior for \5(p\. The \5<p\ trajectory which we get by the 
analytic treatment, as given in Table 1, should be interpreted as an "attrac- 
tor" for solutions. Indeed, if inflation is powerful (large e-fold number), the 
time at which the perturbation is generated becomes unimportant and \8ip\ 
will rapidly approach the behavior found analytically. For the cases where 
inflation is not so powerful (see Fig. 1), the shape of the fluctuations of \S(p\ 
will depend more on the details of the background solution from which they 
were produced. That is exactly what we find in Fig. 3. After inflation the 
long wavelength perturbations of \8<p\ oscillate. 

For an initial fluctuation 5(p, which lies inside the Hubble radius, the 
exponential damping is so effective, that at the end of inflation S(p practically 
vanishes, as can be seen in Fig. 5. This is valid whatever the values of A and 
m are. 

We notice that, to check numerically the analytic behavior for the situa- 
tion where S ~ S max , is more difficult, since one has to deal with numerical 
instabilities. 
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6 Concluding Remarks 



In this paper we studied the perturbation of the coupled Einstein-Klein- 
Gordon equations for a complex scalar field in different regimes. We derived 
analytically the time evolution of long and short wavelength perturbations 
during inflation. Moreover, we established that in order to have a long period 
of inflation the bosonic charge must be close to zero. After inflation the 
perturbation of the metric remains constant during the oscillatory phase. 

If the complex scalar field has a large bosonic charge or the values of the 
parameters entering in the potential are such that the slow-roll approximation 
does not apply, then it turns out that inflation can not have taken place 
or lasted only for a short period. Hence, the background passes directly 
in the oscillatory phase. We showed that in the massive case, for H ~ or 
H ~ H max , the perturbations of the scalar field decay or, at best, stay constant 
and thus no gravitational seeds, like boson stars, can form during this epoch. 
Constraints on the parameters of the scalar field potential coming from the 
measured CMB anisotropics can only be put once a thorough study of the 
quantized perturbations will be done. However, since inflation is essentially 
driven by one component of the complex field, we do not expect significant 
differences to occur on the constraints and thus on the physical implications 



as compared to the real scalar field case ||10|| . 

We also derived the Jeans wavenumber for the Newtonian regime starting 
from the general relativistic equations. 
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Caption 



Table 1: 

Summary of the analytic solutions during the inflationary and the oscillatory 
phases for the metric and the complex perturbations. 

Figure 1: 

h = aH as function of the cosmological time t for three different values of 

AJVf 2 

A = —T- and for the case m = 0, A = 1. A is the coupling constant of the 
quartic self-interaction term, m is the mass of the scalar field and M p is the 
Planck mass. For the massive cases h (t) is in units of m (m _1 ) and for the 
massless case h (t) is in units of M p (M~ v ). 

Figure 2: 

Bardeen potential ip for the zero mode (k = 0) as a function of the cosmo- 
logical time t for different values of m and A. The units of t are the same as 
in Fig. 1. 

Figure 3: 

Norm of \5<p\ for the zero mode as a function of the cosmological time t for 
different values of m and A. For the massive cases (massless case) \8ip\ is in 
units of m (M p ). For the units of t see Fig. I. 

Figure 4: 

Same as for Fig. 2 but with k = 60000 
Figure 5: 

Norm of \S(p\ for k = 60000 and A = as a function of the cosmological time 
t. For the units of \5ip\ and t see Fig. 3 and Fig. 1. 
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During Inflation 




A = 


A ^ 


k 2 > g/g 

(short wavelength) 


ip: oscillates 
Sip: exponentially 
damped 


■0: oscillates in an expo- 
nentially decaying 

envelope 
5tp: exponentially 

damped 


k 2 < g/g 

(long wavelength) 


^ oc w 
\5<p\<x-\ 
(t -> -oo) 


oc exp (4M pV /A/3t) 

\5ip\ oc exp (2M p ^A/3t) 
(t -> -oo) 



During the Oscillatory Phase 




m/0 and E/E max « 1 


k 2 < |Ai - /i 2 | < p'/p 
and 

\h - h 2 \ < A; 2 < p'/p 
(long wavelength) 


^ oc 1/a oc l/t 2/3 

I'M oc | 
(t -> +oo) 


A; 2 > p'/p 

(short wavelength) 


ip: 2 terms, the first oscillates in a | en- 
velope, whereas the second is propor- 
tional to l/t 2/3 for t — > +oo 
2 terms, the first oscillates and the sec- 

/i \ 2 /3 

ond oscillates in a (-H envelope, for 
t — > +oo. 



Table 1 
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